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Abstract. - We analyse the universal properties of nonequilibrium steady states of driven 
Magnetohydrodynamic (MHD) turbulence in three dimensions (3d). We elucidate the depen- 
dence of various phenomenologically important dimensionless constants on the symmetries of 
the two-point correlation functions. We, for the first time, also suggest the intriguing possibility 
of multiscaling universality class varying continuously with certain dimensionless parameters. 
The experimental and theoretical implications of our results are discussed. 



In the vicinity of a critical point, equilibrium systems show universal scaling properties 
for thermodynamic functions and correlation functions: These are characterised by universal 
scaling exponents that depend on the spatial dimension d and the symmetry of the order 
parameter (e.g., Ising, XY etc.) |T], but not on the parameters that specify the details of the 
Hamiltonian. A notable exception is the class of two-dimensional models, such as the XY 
model in d = 2, in which continuously varying scaling exponents are found PJ. Dynamic scal- 
ing exponents, that characterise the behaviour of time-dependent correlation functions also 
show universality A different situation arises in driven, dissipative, nonequilibrium systems 
with nonequilibrium statistical steady states (NESS). We show here that driven homogeneous 
and isotropic MHD turbulence in three dimensions (3d) is a good natural candidate for a 
system with an NESS whose universal properties vary continuously with the degree of cross- 
correlations between the velocity and magnetic fields. Our results illustrate the sensitivities 
of the statistical properties of nonequilibrium steady states on the parameters of the model, 
a situation which does not arise in equilibrium systems. 

Turbulence in plasmas is described by the equations of MHD 4 in 3d for the coupled 
evolution of the velocity field v and the magnetic field b. Solar-wind data and numerical 
simulations of the 3dMHD equations as well as shell models [B] show multiscaling corrections 
to the simple Kolmogorov (K41) scaling [7] for the structure functions of v and b: S%(r) = 
([|cii(x + r) — ai(x)| p ]) ~ r^p , a = v,b for r in the inertial range between the forcing scale L 
and the dissipation scale rjd (i-e., L ^> r ^> In observations and experiments magnetic 
Prandtl number P m (=the ratio of the magnetic to kinetic viscosity) and e = -^p- (the ratio 

(*) email:basu@hmi.dc 
© EDP Sciences 



2 



EUROPHYSICS LETTERS 



of the magnetic- to the kinetic- energy) are known to have a wide range of values: e <C 1 => 
kinetic regime of 3eflVlHD (as on the solar surface), e ~ O(l) =>■ equipartition regime (e.g., 
solar wind in its rest frame) and e>l ^ magnetic regime (e.g., fusion plasmas) 8 . Similarly 
P m can be high (> 1), e.g., in galactic and proto galactic plasmas (Hj. Low values of P m are 
found in liquid metals |H|. In numerical simulations, the structure functions of v and b in 3d 
exhibit multiscaling with exponents different from their fluid analogues |fil lll| . In particular 
Ref. elucidated different multiscaling exponents for v and b (i.e., £p ^ Cj). 

In this letter we develop a one-loop, self-consistent theory for the randomly forced 3dMHD 
equations and use it to calculate several interesting results concerning its NESS. We show, 
in particular, that (i) the variation of e and renormalised (effective) P m are connected to 
each other (ii) P m , e, Kolmogorov's constants K v and K b (see below) vary with the degree 
of crosscorrelations between v and b. We calculate the intermittency exponents (these in 
a log-normal model provide approximate estimates of the deviations from the K41 scaling) 
which also vary continuously with the degree of crosscorrelations between v and b. This 
suggests an interesting possibility of continuously varying multiscaling universality classes for 
3efiVlHD. We show that these variations are linked with the ratios of the different two-point 
correlation functions. Our observations on multiscaling of 3dMHD based on our self-consistent 
calculations on a log-normal model allows for different multiscaling exponents for the velocity 
and magnetic fields as reported in Ref. |6]. 

The 3dMHD equations consist of the Navier-Stokes equation for the velocity field v sup- 
plemented by the Lorentz force due to the magnetic field b: 

<9v „ Vp (Vxb)xb o 

— +v.Vv = + ± —!■ + ^ V 2 v + f 1 

at p inp 

with V.v = (incompressibility) and Ampere's law for a conducting fluid 

9b 

— +v.Vb = b.Vv + /i V 2 b + g. (2) 

In and J5J, p is the fluid density, p is the pressure, p is the "magnetic viscosity", 
arising from the nonzero resistivity of the plasma, v is the kinematic viscosity and f and 
g are forcing functions. Under rescaling of space and time v(x, t) and b(x,i) scale as 
Vi(lx,l z t) — > l Xl Vi, bi(lyi,b z t) — » l X2 bi(x,t). Since v and b are polar and axial vectors respec- 
tively, crosscorrelation tensor (t>j(k, t)bj(— k, 0)) is odd and imaginary in wavevector k |12j . 

We employ a self-consistent mode coupling scheme (SCMC). The response and the correla- 
tion functions of the fields v and b are defined by Vi = G v fi, bi = Gbgi, and correlators = 
{Vi{k,t)vj(-k,0)), C| = (bi(k,t)bj(-k,0)) with the crosscorrelator = (^(k, t)bj (-k, 0)). 
In the scaling limit, in terms of the dynamic exponent z and the roughness exponents xi an d 
X2, G" 1 = iu — T, v {k,w), G^ 1 = iu — £(,(fc,cj) where E„ = k~ z ri v (uj/k z ), Sfc = k~ z r] b (uj/k z ) 
are the self-energies whereas the correlators are given by Q) = Pijk~ d ~ 2xi ~ z rii(uj/k z ), Cfj = 
Pijk~ d ~ 2x2 ~ z rj2{ijj/k z ). rj v (uj/k z ), ribioj/k 2 ), r]i(u;/k z ), T]2(co/k z ) are scaling functions and Pij 
is the transverse projection operator. In systems out of equilibrium, there is no particular 
relation between the noise variance and the dissipation coefficient, unlike their equilibrium 
counterparts where such a relation exists due to the Fluctuation-Dissipation-Theorem. It is 
also well-known that the statistical properties of the nonequilibrium systems depend strongly 
on the noise variances. We, here in particular, take the forcing terms f, g to be Gaussian 
noises with zero mean and covariances proportional to k~ d in d-dimension |13) . This ensures 
that, in absence of a mean magnetic field, the energy spectra are K41-like |fil7| . In this model, 
the advective nonlinearities and the noise variances do not renormalize. This means under 
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simultaneous rescaling of space and time, the nonlinear vertices and the noise variances are 
affected only by their naive scaling dimensions. At the fixed point, these conditions leads to 
strong dynamical scaling, i.e., same dynamic exponent for both the fields and immediately 
yield z — 2/3, \i — X2 — X = 1/3 corresponding to K41 scaling ^JE]- Hence the ratios 
of various correlations are dimensionless numbers. We assume Lorentzian line shapes for the 
self-energies and the correlation functions. The fact that the values of the scaling exponents 
z and x are dimension independent is because the noise correlations change as a function of 
dimensionality in such a way as to render the exponents dimension independent. 

We start with the zero-frequency forms for the self-energies (or the relaxation rates) 
E v (k,co = 0) = vk\ E b (k,0J = 0) = fik z , and the correlations C£(fc,w = 0) = ^i^-Ar 2 * - '*-*, C^(k, to = 
0) = ^jjr Pijk~ 2x ~ z ~ d , in d-dimension. The crosscorrelation tensor Cfj in general has both 
symmetric C?- and antisymmetric C?- parts: C?- = C?-+C£; Cf J (k,uj = 0) = 2iD M p a k -2 X -z-d^ 
ind-dimensionandwithZ)(k) = -D(-k); D(k)D(k) = (D) 2 and C? (k, uj = 0) = ^§-e lJp k p k- 2 ^- z - 4 
(this form is for 3d only). The mean crosshelicity X)k( v (k)-b( — k)) an( i mean electromotive 
force X)k( v (' < ) x M — -K)) are proportional to D and D respectively. Here e« p is the totally an- 
tisymmetric tensor in 3d. Previous renormalisation group/self-consistent approach to 3dMHD 
did not take any crosscorrelation function ^3] fully into consideration. The noises in stochsti- 
cally driven models, such as Eqs. Q and (J2J, represent coarse-grained effects of the degrees 
of freedom whose time-scales of decay are faster than the coarse-grained variables (v and b) 
used. Since both the noises in Eqs. (QJ and 101 have same microscopic origin there is no apri- 
ori reason for their crosscorrelations to vanish always. In SCMC approach vertex corrections 
are neglected. Lack of vertex renormalisations in the zero wavevector limit in 3dMHD allows 
SCMC to yield exact relations between the scaling exponents z and x, as in the the noisy 
Burgers/Kardar-Parisi-Zhang equation [T^. In the context of the noisy Burgers equation in 
I+f dimensions, Frey et al have shown, by using nonrenormalisation of the advective nonlin- 
earities and a second-order perturbation theory, that the effects of the vertex corrections on 
the correlation functions are very small |16j . Similar conclusions for 3dMHD will presumably 
follow from the vertex nonrenormalisations in 3dMHD, though a rigorous calculation is still 
lacking. 

One may note that the symmetries of the symmetric and the antisymmetric parts of the 
crosscorelations are different under the exchange of the cartesian indices. Hence, the presence 
of one does not lead to the generation of the other by non-linear interactions. Thus it is 
possible to treat them separately without any loss of generality. It is straight forward to 
extend the various expressions derived in the paper for the case when both the symmetric and 
the anti-symmetric parts are present. 

Symmetric Cross-correlations.- We obtain the one-loop self-consistent integral equations 
for the self-energies and correlation functions. The one- loop diagrammatic corrections to the 
crosscorrelation function vanish as it is odd in momentum (see, e.g., Q7]). We first consider 



the case when D = 0, i.e., no antisymmetric part of the crosscorrelations. By matching at 
u> = 0, for the self-energies we obtain (by using y = 1/3, z = 2/3) in any dimension d 



v = 



3Di S d d-1 
~4v 2 (2it) d (d + 2) 



+ 



3£> 2 S d d 2 + d-4 
V (2-K) d d(d + 2) ' 



(3) 



3£>i S d 




■3D 2 S d 




(4) 



V = 



2v(v + (i) {2n) d 



2\x(y + fi) (27r) d 
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Similarly, by demanding consistency in the amplitudes of the correlations we get 



r d 2 


Dl 4D 2 




1 


2 , 2 








d ~r d(d+2) 




4DiD 2 1 16D 


2 


(1-1) 




vfj.(u+fj.) 3 {y-\ 





(5) 



In terms of the renormalised magnetic Prandtl number P m = ji/v and V = D2/ D\ we can 
write Eqs.©, Q and © as 



rf-1 
2(d+2) 



r 

2P 2 



2 +d-4 



d(d+2) 



1 + P, 



(1-^ + 7-1^(1-!)' 



(6) 



and T 



r 

P':i 



4/3 



P m (l + P,r 



1 _ 2 j 2_ 

± d ' d(d+2) 



4r-i 



16/3 



p m (i+p m ) T (i+P„) 3 r 2 



(1-1) 



(7) 



where (3 = {jy^j 2 ■ Equations ||BJ and J7J) explicitly demonstrate that P m and the magnetic- 
to kinetic- energy ratio e = ■£- depend on each other. When = 0, i.e., zero crosshelicity, 
Eqs.® and Q yield P m = 0.67 and T = 0.63 and hence energy-ratio e = T/P m w 0.94 in 3d. 
Thus for zero crosshelicity P rn and e are fixed numbers. But with finite crosshelicity, i.e., a 
finite (3 in Eqs.© and J7J, P m and T and hence e depend on /3 in any dimension d. For /? > 
few representative values of P m and Y in 3d are given in Table 1. 

Anti- Symmetric Cross-Correlations. - For the effects of the antisymmetric part (under the 
exchange of i and j, i,j are Cartesian indices) of the crosscorrelation function in 3d on the 
long wavelength properties we again use self-consistent methods. Equations and J3J remain 
unaltered. However, Eq.JSJ) is changed to 



Di 
Do 



' n 2 n 2 
—4- + 2 
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■ 4D 2 
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16D 2 








3(v+fj.) 3 





(8) 



Again by defining /3 



' IL 



) 2 , we obtain self-consistently 



4/3 
3(1+P, r 



4r- 



16/3 



3P m (l + P m ) 3(1 + P m ) a 



(9) 



Thus Eqs.© and @ together provide self-consistent relations between P m and T and hence 
between P m and e when there is a finite antisymmetric crosscorrelation. A few representative 
values of P m and T for /3 > are given in Table 1. 

Kolmogorov's constants.- According to the Kolmogorov's hypothesis for fluid turbulence 
in the inertial range energy spectrum E(k) = K e 2 / 3 k~ 5 / 3 , where K , a universal constant, is 
the Kolmogorov's constant and e is the energy dissipation rate per unit mass. Various calcu- 
lations, based on different techniques by different groups 18 19, 13,20 show that K a ~ 1.5 in 
three dimensions. For MHD, by using Novikov's theorem we connect the dissipation of the 
total energy with the sum of the velocity and magnetic field correlation-amplitudes. We can 
then relate the dissipation of the reduced energy E^ = E v —Ei> with the correlation- amplitudes. 
Thus, in the lowest order of the perturbation theory we calculate the coefficients of the velocity 



Abhik Basu : Statistical properties of magnetohydrodynamic turbulence 



5 



and magnetic field correlations in terms of dissipation rates: E v (k) = K v <?J Z k 5 / 3 ; E b {k) = 
K b e 2 b /3 k~ 5 / 3 . K v and K b are the Kolmogorov's constants for MHD, and e v and e& are the 
mean dissipation rates of the kinetic energy and magnetic energy respectively. In absence 
of a mean magnetic field energy spectra scale as fc~ 5 / 3 in the inertial range [Hj- In Ref. [22] 
Kolmogorov's constants for MHD turbulence (for the Elsasser variables) have been calcu- 
lated incorrectly for various energy-ratio e due to the absence of crosscorrelation functions of 
appropriate structure. We work out the Kolmogorov's constants with finite crosscorreletion 

S 3 ■ ltfl i 4 T 



functions. From Eq.JBJ) in dimension d — 3, we have v — f 0,6 ^ 1 r§^s \ * (1 + f-pr) ' ■ In our no- 



tations, («j(k, ui)vj(— k, — w)) = 2 ^-y^ik iz wn i cn m 3d gives (in the inertial range) for kinetic 

energy spectrum E v {k) = 1.186[2£>i(l + .62{1 + - Pm & Pm) ^ (2*V (i+i^J^i/s ■ To 

connect dissipation rate e v with noise correlation amplitude D\ we use Novikov's theorem |21) 
which gives 2D\ T^rp = £v Substituting this in E v (k) above, we obtain Kolmogorov's constant 

[1 + ,62{1 + g- - p (1 4 ° J ] 2 / 3 
= 1.186- 5= — 1 Pm(1 / +P '" ) , (10) 

where a = (3, —(3 for finite symmetric and antisymmetric crosscorrelations respectively. A 
similar calculation yields K b — K v . Our results show, unsurprisingly, that K v and K b are not 
fixed numbers, rather they depend upon other dimensionless numbers like magnetic Prandtl 
number P m and energy-ratio e. For T = 0, i.e., for no magnetic fields we get the result 
K v = 1.63 for pure fluid turbulence which is well- within the range of accepted values ^9 23 . 

Possibilities of variable multifractality . — Experiments and numerical sumulations [Hllll| 
find nonlinear multiscaling corrections to the K41 prediction of £° = p/3 for the structure 
functions in the inertial range. Until the date, no controlled perturbative calculation for 

is available. To account for multiscaling in fluid turbulence, however, Obukhov |23] and 
Kolmogorov 1^ assumed a log- normal distribution for dissipation e to arrive at Sp(r) = 

(\Av\p) = C p e p/3 rP/ 3 (£) (,5/2)l,(p - 3) j where z is the mean value of e and ( c ( x + r )e(x)) cx 

((At>) 6 /r 2 ) ~ (L/r) 9 . For small 5, 5 ~ 95. A standard calculation on the randomly stirred 
model yields intermittency exponent 5 — 0.2 |26| where 5 = 95. whereas the best possible 
estimates from experiments is 0.23 |26| . In MHD dissipations of kinetic- (e„) and magnetic- 
(eft) energies fluctuate in space and time. Consequently we define two intermittency exponents 

S v , 5b for the kinetic [e^ = \v^2i a p (§fj + ZST") L anc ^ magnetic energy dissipations [e;, = 

|/i^ a(5 (f^j + \ ] to explain the deviations from the K41 scaling. We calculate the 

exponents below by using the Eqs.([[J and following closely Ref. |2H|. We work with the 
self-consistent forms for the self-energies and correlation functions given above along with the 
consistency relations for the amplitude-ratios T,/3,$ and P m . 

Following Ref. (22|, we find the dissipation correlation functions in 3d to be (e„(x + r)e„(x)) - 
12.4e 2 a 2 i^ 2 ln^, (e b (x + r)e 6 (x)> ~ 12.4e 2 a 2 r 2 i^ 2 In ^. a v is defined by the relation v = 

a v el /3 . By using Eqs.©, Q, @, and © we find, in 3d, 

a v = 0.4[1.62 + 0.62(-J - ^ )?/*[! + 1.33-L] 1 / 3 , (11) 

where a — (3 and —p for finite symmetric and antisymmetric crosscorrelations respectively. 
For the pure fluid case a v ~ 0.5 ^3] and is universal. For MHD, however we see that a v 
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varies with j3 and (3, similar to K v . We then obtain fluid intermittency exponent 




r 2 4a 



) 



-1 4/3 



s v = 



p3 p„h + p\ 



(12) 



where, a is the same as before. A similar calculation obtains magnetic intermittency expo- 
nent <5t, = T 2 S V . As expected, both S v and 5^ vary with P m and T (or e). In general, the 
intermittency corrections to the simple K41 scaling of the v and b structure functions are 
unequal, as has been seen in the direct numerical simulations (DNS) and shell-model stud- 
ies UJ. The presence of multiple intermittency exponents makes it difficult to directly extend 
the fluid log-normal intermittency model to 3dMHD. However, the continuous variations of 
the intermittency exponents with (3 and (3 strongly suggest continuously varying multiscaling 
universality classes for 3dMHD. 

In our calculations, f3 and parametrise the internal symmetries of the system. They 
represent the relative strengths of the symmetric and anti-symmetric parts of the crosscorre- 
lation function. Therefore, (3 and (3 charactrerise the symmetries of the two-point correlation 
functions under inversion of the parity and exchange of the cartesian indices. Results from our 
log-normal model-type calculations open up an intriguing possibility of multiscaling properties 
varying continuously with symmetries parametrised by the dimensionless numbers described 
above. We believe that our results (despite the limitations of mode coupling theories - see, 
e.g., |13H27| 1 will provide a valuable insight to the understanding of the difference in the 
multiscaling properties of the velocity and magnetic fields |B] and may also explain the some- 
what different numerical estimates of the multiscaling exponents by different groups |6I11| . We 
believe our results will stimulate further numerical as well as experimental studies for the mea- 
surements of the dimensionless parameters introduced above. To check our results we suggest 
MHD experiments on liquid sodium systems with active grids and random passage of current. 
Various (3 and (3 can be achieved in experiments/simulations by controlling the appropriate 
Grashof numbers |28| , constructed out of the various noise correlators in stochastically driven 
MHD and by measuring the cross helicity, the electromotive force and the kinetic energy in 
the steady state. Renormalised P m can be calculated by measuring time-dependent correla- 
tions. In closing we note that many natural systems, e.g., solar wind, plasmas in tokamac 
have strong mean magnetic fields making them anisotropic. Such effects can also be included 
in our scheme of calculations by appropriately modifying the form of the self-energies. 

From a broader point of view, our results unveil qualitatively new properties in the steady- 
states of nonequilibrium driven, difusive systems. In contrast to the more well-known equilib- 
rium systems, we illustrate the richness of nonequilibrium phenomenologies by using 3dMHD 
as an example. Recently, Drossel et al (2H| , in a set of coupled Langevin equations describing 
the interplay between phase ordering dynamics in the bulk and roughening dynamics of the 
interface of binary films, find a similar continuous variation of the dynamical exponent with 
the coupling strength of the bulk and surface fields. Recently, continuously varying univer- 
sal properties in the context of a simple coupled Burgers-like model has been discussed in 
Ref. I3()| . Here, for the first time, we illustrate a natural realisation of continuously varying 
multiscaling universality class in 3dMHD. We believe our studies will provide further stimuli 
to the studies of these new aspects of driven nonequilibrium systems. 

The author thanks the Alexander von Humboldt Stiftung, Germany for financial support 
and R. Pandit, E. Frey and J. Santos for their critical comments on the manuscript. 



Abhik Basu : Statistical properties of magnetohydrodynamic turbulence 



7 



Table I - Representative values of P m , F, e, K v , ct v and S v as functions of (3 and (5. 
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S v 
















1.6 


0.47 


0.2 
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1.76 


1.76 
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0.15 


3.5 


1.6 
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0.67 


0.63 
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